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Abstract 

Genetic diversity was assessed in samples of cultured Atlantic salmon, Salmo salar L, obtained from facilities in Chile 
between 2005 and 2010, a period of time during which the infectious pathogens Infectious Salmon Anemia (ISA) virus, 
Callgus rogercresseyi (sea lice), and Plscirickettsia salmonis (salmon ricl<ettsial syndrome) were common. Two panels of 
microsatellite markers were utilized: one with microsatellites with no known gene associations (neutral) and one featuring 
microsatellites linked to putative immune-related genes (immune-related). Allelic richness and gene diversity across samples 
were significantly greater in neutral loci as compared to immune-related loci. Both diversity measures were homogeneous 
among samples for immune-related loci and heterogeneous among samples for neutral loci. Immune-related loci were 
identified as Fst outliers in pairwise comparisons of samples at a 10-fold higher frequency than neutral loci. These results 
indicate that neutral and immune-related portions of the Atlantic salmon genome may have differed in response to the 
gauntlet of pathogens and that monitoring of specific, well characterized immune-related loci as well as neutral loci in 
cultured species could be useful when disease control and prevention is a goal. 
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introduction 

Genetic assessment has become a valuable component of many 
aquaculture programs. Keeping track of parentage, marker- 
assisted selection, and monitoring standing genetic variability of 
brood stocks are among the more important applications [1]. 
Selectively neutral, single-copy markers such as microsatellites 
have been used widely because their design and implementation 
are fairly straightforward, allowing comparisons between wild and 
captive populations [2-4]. In addition, measures of variability 
obtained with large panels of microsatellites often are taken as a 
proxy for population-wide fitness [5-6]. The strength of the 
correlation between neutral genetic variation and aspects of 
quantitative genetic variation associated with fitness, however, is 
controversial [7-8] . In settings with differential selection pressures, 
levels of variation may differ between parts of the genome which 
are selectively neutral and those associated with adaptive fitness 
[9]. This could be a concern in aquaculture facilities where 
unintended selective regimes caused by pathogens and parasites or 
unanticipated mass mortality may alter adaptive parts of the 
genome yet leave neutral parts unaffected [10]. 

Disease management in the culture of Atlantic salmon, Salmo 
salar L., has become a particularly important issue worldwide as a 
number of potentially devastating pathogens, sometimes in 
concert, have impacted the global industry. Examples include 



infectious pancreatic necrosis virus or IPNV [11], salmonid 
alphavirus or SAV [12], infectious salmon anemia virus or ISA 
[13-14], and a variety of parasites [15]. A recent ISA epizootic in 
Chile, which began in 2007, was particularly devastating to the 
Chilean salmon industry, which by 2005 had grown to be the 
second largest producer of Adantic salmon in the world [16]. 
While initial mortalities attributable to the virus were documented 
in 2007, the virus likely was present as many as eleven years prior 
to its detection [17]. To comphcate matters further, sea lice [Caligus 
rogercresseyi) were prevalent across the Chilean aquaculture industry 
[18] and may have increased susceptibility of fish to ISA and 
facilitated spread of the virus [19-20]. In addition, salmon 
rickettsial syndrome, caused by the intracellular bacteria Phcirick- 
ettsia salmonis, had been an issue in Chilean aquaculture facilities 
since the late 1980s [21]. The net result of these pathogens across 
the industry in Chile was a reduction in production of 
approximately two-thirds due to direct mortality and/or culling 
[22]. Because of the exogenous origin of the ISA virus [17], part of 
the recovery plan banned the importation of salmon ova [23]; 
consequendy, present-day brood stocks in Chile are an amalgam 
of shared individuals that survived the epizootic. By 2010, when 
the industry began to recover, there was a loss of ~30% of export 
value and a ~50% reduction of the associated work force [23]. 

We assessed genetic diversity within and among samples from 
seven different farms in Chile. The samples were part of a tissue 
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Table 1. Panels of neutral loci and immune-related loci used to screen samples of Atlantic salmon, Salmo salar, from the salmon 
industry in Chile. 



Locus Type 


Motif" 


Range** 


Ta' 


Label" 


Primers* 


Reference* 


Neutral 


SSsp 1605 


(GATA),, 


214-252 


58 


Fam 


F: CGCAATGGAAGTCAGTGGACTGG 
R: CTGATTTAGC 1 1 1 1 1 AGTGCCCAATGC 


[50] 


SSsp2201 


(GATA)34 


251-335 


58 


Fam 


F: TTTAGATGGTGGGATACTGGGAGGC 
R: CGGGAGCCCCATAACCCTACTAATAAC 


[50] 


SSsp2210 


(GTTA),, 


112-168 


58 


Fam 


F: AAGTATTCATGCACACACATTCACTGC 
R: CAAGACCC 1 1 1 1 1 CCAATGGGATTC 


[50] 


SSsp2213 


(GTTA)22 


151-215 


58 


Hex 


F: ATGTGGAGGTCAACTAACCAGCGTG 
R: CATCAATCACAGAGTGAGGCACTCG 


[50] 


SS5p2215 


(GTTA),4 


107-175 


58 


Ned 


F: ACTAGCCAGGTGTCCTGCCGGTC 
R: AGGGTCAGTCAGTCACACCATGCAC 


[50] 


SS5p2216 


(GTTA,25 


201-257 


58 


Ned 


F: GGCCCAGACAGATAAACAAACACGC 
R: GCCAACAGCAGCATCTACACCCAG 


[50] 


SSspG7 


(GTTA),8 


121-215 


58 


Hex 


F: CTTGGTCCCGTTCTTACGACAACC 
R: TGCACGCTGCTTGGTCCTTG 


[50] 


SSa85 


(GT),4 


106-186 


58 


Ned 


F: AGGTGGGTCCTCCAAGCTAC 
R: ACCCGCTCCTCACTTAATC 


[51] 


SSa 197 


(GT)5C(TG)4TC(TG)3A(GTGA),5 


162-260 


58 


Fam 


F: GGGTTGAGTAGGGAGGCTTG 
R: TGGCAGGGATTTGACATAAC 


[51] 


SSal 71 


{TGTA),4(TG)7 


198-252 


58 


Fam 


F: TTATTATCCAAAGGGGTCAAAA 
R: GAGGTCGCTGGGGTTTACTAT 


[51] 


SSa202 


{CA)3(CTCA),7 


237-371 


58 


Ned 


F: CTTGGAATATCTAGAATATGGC 
R: TTCATGTGTTAATGTTGCGTG 


[51] 


SSal 4 


(TC),oN,5(TC)3N2(AC),2(TC)3N5(CA)4 


161-165 


NA 


Tail 


F: CCTTTTGACAGATTTAGGATTTC 
R: CAAACCAAACATACCTAAAGCC 


[52] 


SSOSL85 


(GT)22 


201-241 


NA 


Tail 


F: TGTGGATTTTTGTATTATGTTA 
R: ATACATTTCCTCCTCATTCAGT 


[53] 


SS0SL311 


(TG)38 


139-190 


NA 


Tail 


F: TAGATAATGGAGGAACTGCATTCT 
R: CATGCTTCATAAGAAAAAGATTGT 


[53] 


SSOSL438 


{AC)26AT(AC)8 


132-174 


NA 


Tail 


F: GACAACACACAACCAAGGCAC 
R: TTATGCTAGGTCTTTATGCATTGT 


[53] 


Locus Type 


MotiP 


Range'' 


Ta' 


Label" 


Primers' 


Reference' 


Immune-related 


SsalROOlTKU 


(AG), 2 


113-115 


NA 


Tail 


F: AAGAGCGAGAGAGAAGGATGG 
R: GTTTCACAGAATCAACAGTCAGCAA 


[47] 


SsalR002TKU 


(AG),o 


241-249 


NA 


Tail 


F: GGGTACAAGCAGGGGTCTTA 

R: GTTTAAGAGTGGACCGACAACAAT 


[47] 


SsalROOSTKU 


(GT),3 


146-163 


NA 


Tail 


F: TGTTCTGCAGGTCAGAAGTGA 

R: GTTTGAGTGGGAGGAAGGGGAGTA 


[47] 


SSOIR004TKU 


(AG),o 


313-321 


NA 


Tail 


F: AGCTATTTCCAAGGCGTTCA 

R: GTTTCACCACTCAGGAGAGCATGA 


[47] 


SsalROOSTKU 


(CTT)7 


371-379 


NA 


Tail 


F: CGACGACTTTTTCATCTGTCTT 

R: GTTTGGACAACACATCTCATTCCAA 


[47] 


SsalROOTTKU 


(AG)9 


183-191 


NA 


Tail 


F: GCAATGCTGCCATCTAGTGA 

R: GTTTCAAGGAAAGCCTACAAAAAGC 


[47] 


SsalROlOTKU 


(AG),o 


192-198 


NA 


Tail 


F: CAACGACACCATACCAACCA 
R: GTTTAACCCCTTCCAAGTTCCATC 


[47] 


SsalROllTKU 


(ATC)8 


364-376 


NA 


Tail 


F: CCAGCCAACTACGACAACTG 
R: GTTTGTGGTTATTTTTGGGGTGA 


[47] 


SsfllR012TKU 


(AATC),4 


162-222 


NA 


Tail 


F: GAGTCCCCTTTGGCCTCTC 

R: GTTTAAACACAGTAAGCCCATCTATTG 


[47] 


SsalROMTKU 


(ATOs 


395-407 


NA 


Tail 


F: CTGAGGTGGTGGCACAGC 

R: GTTTATTGTTTGGTTCTTACAGCAGGA 


[47] 


S5fllR016TKU 


(AT), 


317-327 


NA 


Tail 


F: CCAAAAATGTCCCATTCACC 
R: GTTTGTGTGCCACTCAGAATTG 


[47] 
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Table 1. Cont. 







Locus Typ6 


MotiP 




' A 


Label** 


Primers^ 




fiG934281 


(TCTG),4 


215-277 


NA 


Tail 


F: ACTGCTTCTCCCCTGCTACA 

R: GTTTGCGAACCACACATATACCAC 


[46] 


C4769358 


(AC)22 


133-173 


NA 


Tail 


F: TGACGCCATATGCAAAGAGA 
R: GTTTCCTTTGTCTGCAAAACGTGA 


[46] 


UBA 


(CA)17 


330-364± 


NA 


Tail 


F: GGAGAGCTGCCCAGATGACTT 
R: CAATTACCACAAGCCCGCTC 


[44] 


MHCJI 


GTTATTAAAT 


232-292 


NA 


Tail 


F: GATGGCAAAGAGGAAAGTGAG 
R: TTGTTATGCTCTACCTCTGAA 


[45] 



^Motif indicates repeat motif. 

"^Range refers to alleles thus far uncovered (may include the 21 bp 50-tail-sequence}. 

*^Ta is optimized annealing temperature. For loci amplified with the tail protocol, a step-down protocol was used so annealing temperature is not applicable (NA). 
^Label refers to dye {6-Fam, Ned, or Hex) used to label forward primer. Tail refers to loci that were amplified using the tail protocol. 
*^Primer sequences are forward {F and reverse (R). 
'citations describing putative function. 
doi:1 0.1 371 /journal.pone.0099358.t001 



bio-bank maintained at the University de Santiago de Chile and 
procured between 2005 and 2010. Two panels of microsatellite 
markers were utilized: one featuring microsateUites not known to 
be associated with functional genes (neutral loci), and the other 
featuring microsateUites known to be linked to putative functional 
immune genes (immune-related loci). The primary goal of this 
study was to assess whether patterns of genetic diversity within and 
among samples obtained from selectively neutral loci were 
consistent with patterns of diversity obtained from loci which 
potentially could be under selection. 

Materials and Methods 

Tissues from a total of 109 Atlantic salmon were obtained by R. 
Vidal from seven aquaculture facilities in Chile. Samples were 
obtained from live fish by a third-party diagnostic laboratory that 
screens aquaculture farms for infectious agents (Diagnotec S.A.) 
We obtained tissues that had been stored at the University de 
Santiago de ChUe. There is no official animal care protocol in 
Chile. However, salmon aquaculture is an enormous industry 
(second largest salmon aquaculture industry in the world) and 
aquaculture farms are tightly regxilated due to past problems with 
widespread disease. Therefore all farms had necessary permits. 
Because of commercial confidentiality issues, individual companies 
and localities of each facility are not named. Samples from each 
individual source farm were labeled A - G. Samples A - C were 
obtained in 20 1 0 when recovery was underway; whereas samples 
D - G were obtained between 2005 and 2007 during a period of 
high disease-related mortality. Sample sizes ranged from 10 
(sample G) to 30 (sample A). Tissues were taken directly from 
brood stock and preserved in 20% DMSO buffer [24]. DNA was 
extracted with a modified Chelex extraction protocol [25]. After 
centrifugation for two minutes at 13,000 xg, the resulting 
supernatant was used as template for PCR amplification. 

Individuals were genotyped at 30 microsatellite loci; 15 of the 30 
loci assayed were reported in previous studies to behave as 
selectively neutral [26,27,28,29] (hereafter, neutral), while the 
other 15 were reported to be tightly linked to genes relating to 
immune function in salmon (hereafter, immune-related). Informa- 
tion for each microsatellite, including the original reference, is 
given in Table 1 . Each microsatellite was amplified using either a 
fluorescently labeled forward primer or using the 'tail-labeled' 
methodology reported in Karlsson et al. [30] . In both cases, one of 
three fluorescent dyes, 6-Fam, Ned, or Hex (Dye Set D, Applied 



Biosystems), was used. Polymerase chain reaction (PCR) amplifi- 
cation, using fluorescently labeled primers, was carried out in 
10 |J,L reactions containing: Ix reaction buffer, 2 mM MgClj, 
0.25 niM of each dNTP, 5 picomoles of each primer, 0.05 U/^L 
of Taq polymerase, and 50-150 nanomoles of DNA template. 
Reaction conditions consisted of an initial denaturation at 95°C 
for 3 min, followed by 38 cycles of 95°C for 30 sec, 58°C for 
45 sec, and 72°C for 60 sec, and followed by a final extension at 
72°C for 10 min. PCR amplification using the tail-labeled 
methodology was carried out in 10 |J.L reactions containing: Ix 
reaction buffer, 2 niM MgClj, 0.25 mM of each dNTP, 5 
picomoles of the tail and reverse primers, 0.5 picomoles of 
forward primer, 0.05 U/|.tL of Taq polymerase, and 50—150 
nanomoles of DNA template. Reaction conditions consisted of an 
initial denaturation at 95 °C for 15 min, followed by 14 cycles of 
95°C for 30 sec, 59°C for 90 sec, and 72°C for 60 sec and 24 
cycles of 95°C for 30 sec, 52°C for 90 sec, and 72°C for 60 sec, 
followed by a final extension at 60°C for 30 min. Amplicons were 
electrophoresed and visualized on 6% polyacrylamide gels, using 
an ABI 377 Automated Sequencer (AppKed Biosystems). Allele 
calling was performed manually, using Genotyper 2.5 (Perkin 
Elmer) and Genescan 3.1.2 (Applied Biosystems). 

Conformance to expectations of Hardy- Weinberg equilibrium 
(HWE) was evaluated for each microsatellite in each sample, using 
Genepcjp v.4.0 [31-32]; significance was assessed at the 0.05 level, 
using exact tests, with 1,000 batches and 10,000 iterations per 
batch. Sequential Bonferroni adjustment [33] was used to correct 
for multiple testing. Allelic richness (the expected number of alleles 
standardized to the smallest sample size using rarefaction) and 
unbiased gene diversity (expected heterozygosity) were estimated 
for each microsatellite in each sample, using Fst At' v.2.9.3.2 [34— 
35]. WUcoxon signed-rank tests, implemented in Systat 8.0 (SPSS 
Inc.), were used to test for homogeneity in allelic richness and gene 
diversity between marker panels (neutral vs. immune-related). 
Friedman rank tests, as implemented in Systat, were used to test 
homogeneity of allelic richness and gene diversity across samples 
for both neutral and immune-related loci separately. Finally, linear 
regression was used to test for a correlation between diversity 
estimates made with both marker panels (neutral and immune- 
related) across samples. 

Losn'AN [36] was used to screen for F^y outliers (candidate loci 
under selection) by comparing observed mean pairwise values 
at each microsatellite between pairs of samples; mean Fst values 
were corrected for locus-specific gene diversity (expected hetero- 
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Figure 1. Comparisons of mean allelic richness ± standard error for neutral loci (a) and immune loci (b) for samples taken during 
the epizootic (square) and sample taken during recovery (triangle). Locus names (see Table 1) are abbreviated. Figure created using R 
v.2.1.3.1 [57], using dotplots.errors [58]. 
doi:1 0.1 371 /journal.pone.0099358.g001 



PLOS ONE I www.plosone.org 



4 



June 2014 | Volume 9 | Issue 6 | e99358 



Genetic Variation in Chilean Cultured Salmon 



0.700 




1 
1 
1 






0.600 




1 
1 
1 
1 
1 




A 


0.500 




1 
1 
j 
1 






0.400 




1 
1 

! 






0.300 




1 


* 




0.200 
0.100 
0.000 


1 • 


1 1 

1 * A » 1 , 


i 


-0.100 




1 














-16.000 


-8.000 0.000 8.000 
51ogit(P) 


16.000 



Figure 2. Results of outlier tests: ordinate (pairwise Fst values); abscissa (logit transformed probability values). Heavy dashed lines 
indicate significance at a =0.05; light dotted lines indicate significance at a =0.01. Circles are neutral loci; triangles are immune-related loci. 
doi:10.1371/journal.pone.0099358.g002 



zygosity), against a 95% and 99% confidence interval of Fsx^slues 
(corrected for diversity) generated by simulation. Simulations were 
implemented using the stepwise mutational model with 60,000 
steps and an initial simulation of 60,000 steps run to estimate mean 
neutral Fs,. 

Results 

Genotypes at all microsateUites assayed are in Data Set SI. 
Summary statistics for all microsateUites in each sample are 
presented in Table S 1 . Genotypes at all microsateUites conformed 
to the expectations of HWE, foUowing correction for multiple tests. 
Mean aUelic richness assayed per sample (± S.E.) for neutral loci 
ranged from 5.68±0.34 (Sample F) to 8.34±0.66 (Sample B) and 
for immune-related loci ranged from 3.50±0.35 (Sample G) to 
4.54±0.62 (Sample A). Mean gene diversity assayed per sample 
(± S.E.) for neutral loci ranged from 0.746±0.043 (Sample G) to 
0.840±0.044 (Sample B) and from 0.516±0.073 (Sample B) to 
0.602 ±0.058 (Sample A) for immune-related loci. WUcoxon 
signed-rank tests revealed that both allelic richness (^ = 2.366; 
P =0.018) and gene diversity (Z =2.366; P =0.018) in neutral 
loci were significandy greater than in immune-related loci. 
Friedman rank tests revealed different patterns of genetic 
variabUity within marker types. Both aUelic richness (Q_ 
= 41.461, P<0.001) and gene diversity (Q, = 19.143, P =0.004) 



of neutral loci were heterogeneous across samples; whereas the two 
parameters were homogeneous across samples when immune- 
related loci were considered {Q_ =9.696, P =0.138- and Q_ 
= 7.309, P =0.293, respectively, Fig. 1.). Post-hoc, sign-rank 
analysis of variance (Anova), using just neutral loci, revealed 
significant differences in both allelic richness [F = 36.854; P< 
0.001) and gene diversity {F = 17.749, P<0.001) between samples 
taken between 2005-2007 and those taken in 2010 with a higher 
mean diversity for both measures in the samples from 20 1 0 (allelic 
richness -13.89±1.31 vs. 11.84±1.10; gene diversity — 
0.848±0.030 vs. 0.807±0.031, Fig. 1). No significant correlation 
was found across samples between estimates of diversity based on 
the neutral loci and those based on immune-related loci (aUelic 
richness: [15] = 2.471, P = 0.177, gene diversity: i*' [15] =0.155, 
P =0.710). 

A total of 43 (6.8%) outlier values (at ot = 0.05) were detected 
among 630 pairwise comparisons of samples at each of the 30 
microsateUites markers (Fig. 2). Of these, 13 involved neutral loci; 
the remainder (30 comparisons) involved immune-related loci. A 
total of 22 (3.5%) outliers were detected at ot = 0.01 (Fig. 2); of 
these, two involved one of the 1 5 neutral loci assayed, whUe the 
remainder involved nine of the 15 immune-related loci assayed 
(Fig. 2). Of the 20 Fsr outliers involving the nine immune-related 
loci, 10 had values that exceeded neutral expectations, and 



Table 2. Immune microsateUites that were significant positive outliers, at a =0.01, in pairwise comparisons of samples. 





Locus 


Panel 


Coding Genes 


Function 


Citation 


SSalOOlTKU 


Immune 








SSal004TKU 


Immune 


lnterleukin-1 p 


Mediator of inflammatory response 


[54] 


S5a[012TKU 


Immune 


Neutrophil chemotactic factor 


Involved in inflammatory reactions 


[55] 


CA769358 


Immune 


VHSV-induced protein 


Response to hemorrhagic septicemia virus 


[56] 



Coding genes, and their putative function, to which immune-related markers are potentially linked are indicated. 
doi:l 0.1 371 /journal.pone.0099358.t002 
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Figure 3. Allele frequency distributions (bars, left axis) and unbiased gene diversity (triangles, right axis) for two immune loci that 
were significant positive outliers and for two neutral loci. 

doi:1 0.1 371 /journal.pone.0099358.g003 



nine of these 10 involved comparisons between samples taken 
between 2005 and 2007 and samples taken during 2010. In theory, 
Fgx values that significantly exceed neutral expectations may 
reflect directional or diversifying selection [37-38]. The immune- 
related loci linked to genes that may have been under directional 
or diversifying selection were j'&IOOlTKU, 5SaI004TKU, and 
&aI012TKU, and CA769358; three of which have been inferred 
in other species to be involved in response to inflammation and 
infection by hemorrhagic septicemia virus (Table 2). It is also 
worth noting that at the less stringent (o( = 0.05) level, aU but one 
of the 1 7 FsT values that exceeded neutral expectations involved 
one of these same four immune-related loci. 

Discussion 

We examined patterns of variability in selectively neutral 
microsateUites and in microsatellites closely linked to immune- 
related genes within and among samples of cultured Atlantic 
salmon exposed to a number of pathogens. Part of our goal was to 
assess whether neutral genetic diversity was a suitable proxy for 
adaptive genetic diversity. Tests of homogeneity in allelic richness 
and gene diversity between the two genetic marker types across all 
samples revealed significantly greater allelic richness and gene 



diversity in neutral loci. This finding strongly indicates that the 
selectively neutral loci employed were not appropriate proxies for 
the immune-related loci, and moreover, that there may be 
selective constraints on the diversity of immune-related loci, given 
that strong selection against functional allelic variants can lead to 
reduced variation at closely linked genetic markers [39]. In 
addition, because addition or subtraction of repetitive elements 
can alter gene function, there may be selective constraints on 
polymorphism of microsatellites proximal to coding regions [40— 
41]. Consistent with this, the average range of allele size 
(correlated with number of alleles) within loci differed significandy 
between neutral and immune-related loci (/ [14] = — 3.60, P 
= 0.003), with a higher mean of 64.3 bp for neutral loci versus 
23.3 bp for immune-related loci. Levels of variation in neutral loci, 
measured as both allelic richness and gene diversity, also were 
significandy heterogeneous among samples, with greater levels of 
diversity present in samples taken in 2010 (as the industry 
recovered). This is perhaps due to restocking with a mix of 
individuals from other Chilean aquaculture farms which may not 
have had high levels of diseases resulting in crosses of relatively 
outbred individuals. Levels of variation in immune-related loci, 
however, were homogeneous and there was no significant 
correlation across samples in levels of variation (allelic richness 
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and gene diversity) of the two marker types. This, too, indicates 
that neutral loci are not suitable proxies for immune-related loci, 
failing to capture both patterns and levels of immune variation. 

The two genetic marker t)'pes differed substantially at both ot 
= 0.05 and ot =0.01 in the number of -F^r oudiers detected in 
pairwise comparisons of samples, with far fewer pairwise Fgr 
outliers involving neutral loci than involving immune-related loci. 
At both statistical levels, FgY outiiers that exceeded neutral 
expectations primarily involved immune-related loci, and most 
of these involved pairwise comparisons between samples taken 
between 2005 and 2007 and sampk-s takc-n during 2010. In 
general, Fgr values that significantiy exceed neutral expectations 
are hypothesized to reflect directional or diversifying selection [37- 
38], suggesting that allele-frequency changes at immune-related 
loci might have occurred in response to the gauntiet of pathogens 
present in Chilean aquaculture farms in the mid-to-late 2000s and 
that these changes were independent of overall levels of diversity. 
Consistent with this, diversity in outlier loci is higher for two of the 
loci (&flI012TKU and CA769358), lower for another 
(5'6'aIOOlTKU, Fig. 3a), and equivocal for the fourth 
(SsaI004TKU, Fig. 3b) when samples taken between 2005 and 
2007 are compared with samples taken in 2010. By contrast, 
neutral markers revealed a pattern of increased diversity when 
samples taken between 2005 and 2007 are compared with samples 
taken in 2010 (Fig. 3c-d). Because the original purpose in 
obtaining samples at various aquaculture farms was to have a 
tissue database of salmon brood stock in Chile, no tissue samples 
were a\ailable from the same facilities taken at different time 
periods. In addition, confidc'ntiality issues precluded acquisition of 
further samples. Consequently, while we are able to document 
change in the frequency of several alleles at the four, immune- 
related loci (5'6'aIOOlTKU, 5'&I004TKU, and &flI012, and 
CA769358) where Fgr values significantiy exceeded neuttal 
expectations, we cannot unequivocally demonstrate that the 
changes occurred at the same arjuaculture facility. However, the 
putative functions of three of the four principle immune-related 
loci involve responses to inflammation and/or to hemorrhagic 
viral infection, and two of tiie loci (5'&I004TKU and CA769358) 
have been detected as Fgr outiiers and taken as evidence of 
directional selection in wild populations of Atiantic salmon [29] . 
Consistent with this is the finding that only two Fgr outhers 
[SSspl^Qb and SSsp22\G) involved neutral loci, and in both cases 
estimated Fgr values were less than expected under a model of 
selective neutrality. 

While we cannot unequivocally associate the allele changes 
observed at the immune-related loci with the diseases present; the 
immune-related loci identified in this study may represent a good 
starting point for future genomics research. Given that sequencing 
of the Atiantic salmon genome is near completion [42^3] and 
that expressed-sequence-tags (ESTs) used to design PCR primers 
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